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We apply the Coherent Potential Approximation (CPA) to an ex- 
tended Hubbard model to describe disordered superconductors with 
d-wave pairing. We discuss the pair-breaking effect caused by non- 
magnetic disorder in presence of Van Hove singularity. 

1 Introduction 

The original treatments of the influence of magnetic and non-magnetic disor- 
der on superconductors applied to classic BCS superconductors QXj , has been 
reexamined for superconductors with ci-wave symmetry of order parameter 
[0, [| [|, |5|. On the other hand, in high temperature superconductors the 
distance between the chemical potential and Van Hove singularity was found 
to be relatively small. This has lead to formulation of Van Hove scenario 
for high temperature superconductors which says that optimal critical tem- 
perature is reached when chemical potential passes through the Van Hove 
singularity in the density of states ||. However, doping with charge carriers 
does not only change the density in the system but also smear the density 
of states eliminating its singularities and, specially for anisotropic supercon- 
ductors, introduce the electron pair-breaking phenomenon H, || Q. 

2 Superconductivity in a Disordered System 

We start with the single band Hubbard model with an attractive extended 
interaction which is described by the following Hamiltonian [[/]: 
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In the above rii = + is the charge on site labeled i, fi is the chemical 
potential. Disorder is introduced into the problem by allowing the local site 
energy to vary randomly from site to site, c\ a and q CT are the Fermion 
creation and annihilation operators for an electron on site % with spin a, is 
the amplitude for hopping from site j to site i and finally [Ty is the attractive 
interaction (£/y < 0), between electrons on neighbour sites (i ^ j). 
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Figure 1: (a) Projected ci-wave density of states N d (E) in presence of 
disorder for various 5 and the electron density of states N(E) for the clean 
system, (b) — ImSn(i?) for various band fillings 5. 

Here we will assume, for simplicity, that the random site energy £j has a 
uniform distribution £j G [—5/2, 6/2]. Following the usual way we shall apply 
the Coherent Potential Approximation (CPA) to mean that the coherent 
potential S(-E) = H(i,i]E) ||, in a site approximation, is defined by the 
zero value of an averaged t-matrix T(z, i; E). 

The linearized gap equation in Hartree-Fock-Gorkov approximation 0, [/]] : 
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where T c is the critical temperature, ks denotes the Boltzmann constant, 



G U (E) is an averaged electron Green function which defines the weighted 
density of states (DOS) of ci-wave electron states Nd(E): 
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Hii(E) describes the electron self energy in the normal disordered system and 
T]u = 2(cosfca; — cosky). Examples of projected densities of states Nd(E) for 
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disordered 2D system stem are presented in Fig. la while the corresponding 
self energies are plotted in Fig. lb. Equation 1 can be rewritten in terms of 
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Figure 2: (a) N^iui) versus u - full line '1', the fitting curve c/(u + b) for 
c = 0.95 and b = 0.35 - dashed line '2'. (b) T c versus |E |: '1' the standard 
Abrikosov-Gorkov formula, '2' results obtained from Eq. 5, '3' results from 
numerical solving of the gap equation (Eq. 2) for U = —It. 

Matsubara frequencies u n as 

1 = \E\ V ^ii (^") (4) 
2k B T c ^ 2iu n - Tr£(^ n ) 1 1 

In Fig. 2a we show the density N d {iu) = —lmG d (iuj n )/7r versus imaginary 
frequency u, where E = e/ + iuj for Fermi energy e/ chosen at Van Hove 
singularity (e/ = 0, Fig. la). In the standard treatments [4] it was assumed 
to be constant. However in our case, due to the Van Hove singularity, it 
depends strongly on uj. The corresponding projected density N d {iu) can be 
roughly approximate by a simple formula Nd{iuj) = c/(u + b) where c and b 
are constants (Fig. 2a). Approximating also E(zo>) = — 2|£ |sgn(co>) we get 
from Eq. 2 the analytic pair-breaking formula: 

* G) - * + sk) = ,p (I + p ) - * {I + <" + ^f) • (5) 

where p c = |Eo|/(27rT c ) is a pair-breaking parameter and T c o is the critical 
temperature for a clean superconductor. Note that for large b (b — > oo) 
Eq. 4 transforms into the standard Abrikosov-Gorkov equation [p]]. Fig. 
2b presents the comparison between these analytic formulae as well as the 



3 



numerical results obtained from the gap equation (Eq. 2). Both analytic 
and numerical results (Fig. 2b) show that Van Hove singularity make the 
superconducting system more robust in presence the disorder. 

3 Conclusions 

Analyzing the effect of disorder on disordered d-wave superconductor we have 
found the additional influence of the Van Hove singularity. In the presence 
of a weak disorder 5 < t we observe only small change of the projected 
density of states N d (Fig. la) present in the gap equation (Eqs. 2,4) In 
the same time we observe the rapid degradation of T c , which is connected 
with the pair-breaking effect (Eq. 2). Interestingly, Van Hove singularity 
modifies the standard Abrikosov-Gorkov formula, originally obtained for a 
constant density of states, increasing the critical value of |S | which destroys 
the superconducting phase (Fig. 2b). That result is in qualitative agreement 
with the experimental results on Zn substitutions ||. 
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